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Abstract 



The potential of DA<£NE to explore direct radiative decay (j) — * 7r + 7r"7 is 
studied in detail. Predictions of different theoretical models for this decay are 
compared. We find that it should be possible to discriminate between these 
models at DA<I>NE in one year, even assuming a relatively low luminosity C = 
10 31 cm~ 2 sec" 1 . The influence of the decay 4> — > 7r + 7r"7 on the measurement 
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of total cross section a{e + e — > hadrons) by tagging a photon in the reaction 
e + e~ — ► 7r + 7r~7 is also discussed. 

I. INTRODUCTION 

Investigation of CP violation is the most important physical goal of DA$NE, a high 
luminosity e + e~ collider which operates on the 0(1020) resonance. However, thanks to high 
luminosity, there will be a substantial amount of data which may be used to advance our 
knowledge on low energy hadron dynamics and even contribute to precision electroweak 
measurements 

Recently it was suggested [[|, that the annihilation cross section a(e + e~ — > hadrons) 
at energies below the mass of the <fi resonance may be studied at DA$NE using reaction 
e + e~ — > hadrons + 7. By tagging the photon it is possible to determine the pion form factor 
at the momentum transfer below the mass of the meson There are several possibilities 
to improve further the analysis of @ and in this paper we consider the contribution of the 
direct rare decay — > 7r + 7r~7 to the reaction e + e~ — > 7r + 7r~7. Using the terminology of [Q, 
the direct decay contributes to the final state radiation which, for the purpose of the cross 
section measurement, has to be suppressed by an appropriate choice of cuts on the photon 
and pion angles and energies. One of the aims of the present paper is to find out how the 
contribution of the direct decay affects the analysis of Ref. 0. 

Besides that, the rare decay — > 7r + 7r~7 is an interesting process by itself. As one deals 
here with the low energy limit of QCD, the first principles calculations are not possible and 
one has to resort to various models Since the number of models is flourishing, we think 

that the experiments should distinguish between them. In principle, that can be achieved by 
studying the low energy region of the photon spectrum in the reaction e + e~ — > 7r + 7r~7 PJIUf, 
but it is not an easy task. The reason is that the relative phases of the direct decay — > 
7r + 7r - 7 and the pure QED processes (initial (ISR) and final (FSR) state radiation) are not 
predicted by these models. As a consequence, the sign of the interference termQ appears to be 
to a large extent arbitrary. If one assumes that the interference is destructive, the branching 
ratio of the direct decay becomes very small, BR(0 — > 7r + 7r~7) m 4 x 10~ 5 . Under such 
circumstances, a detailed study of the decay — > 7r + 7r~7 will be rather difficult, requiring 
high statistics and a careful control over efficiencies in order to discriminate between different 
models by fitting the photon spectrum. 

In this paper we report on the implementation of the direct decay — ► 7r + 7r~7 into 
the Monte Carlo event generator for pure QED process e + e~ — ► 7r + 7r _ 7 described in ||. 
Our implementation permits to choose between different models for the decay — > -n^-n^^. 
A clear advantage of having a Monte Carlo event generator for these studies is that it 
allows to keep control over efficiencies and resolution of the detector, fine tuning of the 
parameters and also provides for the possibility to generate realistic distributions where the 



Let us note that for symmetric cuts on the pions angles the initial state radiation does not 
contribute to the interference term because of charge parity conservation. 
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reaction e + e — > 7r + 7r 7 is accompanied by radiation of photons collinear to electrons and 
positrons [@. 



II. THE MATRIX ELEMENT FOR vr+vr 7 FINAL STATE 

The matrix element for the direct decay — > 7r + 7r~7 is parameterized as: 
.M = -ieU(Q 2 , Qq)e lj) tM d fia e a , d m = {Qq)g^ a - q^Q a , 



(1) 



where and e are polarizations of the meson and the photon, respectively; q is the 
momentum of the photon and Q is the momentum of the 0. The function f ( f > (Q 2 ,Qq) in 
Eq.(|l]) is the form factor for the direct decay. Its exact form depends on the chosen model. 

Considering production of the meson in e + e~ collision with the center of mass energy 
squared s = Q 2 and its subsequent decay to 7r + 7r~7 final state, we find: 



Q 2 



F^{Q 2 , Qq)v{p2)l^u(pi)d lJia e a , 



(2) 



where the form factor Fa is defined as: 



9<hf<j> 



Q 2 - M 2 + iMsT a 



(3) 



The coupling constant describes the mixing of the photon and the meson and can be 
determined from the decay width of the meson into electron positron pair. Using 



e e 



3 Ml 



(4) 



and T(0 -> e+e") = 1.3246 ■ 10" 6 GeV, M^ = 1.0194 GeV, one obtains = 7.929 • 
1(T 2 GeV 2 . 

Consider now the QED process e~(pi) + e + (p 2 ) — > 7r + (7i"i) + 7r~(7r 2 ) + 7(g). The initial 
state radiation amplitude reads: 



Mi 



-ie*F n {Q\) , 

-V{P2J 



Qi 



7m (Pl - q)la Ia(-P2 + q)j» 



-2p x q 



-Ipiq 



(5) 



where Q\ = i\\ + 7r 2 and n = 7i"i — 7r 2 . 

For the amplitude of the final state radiation we obtain: 



X 



Q 2 

2*2 +g Q , „ «x (-2vrf-g c 



an 



2ir 2 q 2niq 
Then, the differential cross section for e + e~ — ► 7r + 7r~7 can be written as: 



(6) 
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ddtotai = do-QED + da<j ) , 



(7) 



where gIctqed is the contribution considered in 0] 

dor QED ~ \M isr + M fsr \ 2 , (8) 

and 



da a 



| jW/ + 2Re {M* sr M^} + 2Re {M} sr M^}} (9) 



includes the amplitude of the direct decay. One sees that da^ contains different interference 
terms. For this reason one might expect a significant dependence of the — > 7r + 7r _ 7 signal 
on the relative phases of M.$ and A4i sr + M.f sr - We will show below that this is indeed 
the case. We now describe three different models for the direct decay — > 7r + 7r - 7 that are 
implemented in our event generator. 

1. "No structure" model ||. 

In this case the decay — > 7r + 7r~7 occurs through two subsequent transitions — > f Q j ~~ > 
7T + 7r - 7. The form factor in Eq.(|T|) becomes: 

The coupling constants in the above equation can be estimated by using the information 
on the branching ratio Br(0 — > /07) and on the branching ratio Br(/o — > 7t + tt~). One 
obtains 01: 



\94>fay9fo*+*-\ = (144 ± 15) v /Br(0 ^ / 7)- 

Needless to say, that the region of applicability of this model is restricted to relatively soft 
photons, when the fo meson in the intermediate state is not too far off shell. For this reason, 
when implementing this model into the event generator, we have introduced an additional 
exponential damping factor which suppresses the emission rate for high energy photons || : 

f™ str - -> 1.625 x /^exp {-^}, (11) 

with A = 0.3 GeV [|. 

2. K+K- model 

In this model one also has a two step transition, similar to "no structure" model. How- 
ever, the — > /07 decay amplitude is generated dynamically through the loop of charged 
kaons. The form factor in Eq.(|l]) reads: 

rK+K- _ 94>K+K-9f -K+-K-9f K+K- j ( Qj_\ ( , ?] 

^ 2^m%(m%-Ql-im f0 T f0 ) {m* K 'm* K )' 1 } 

The coupling constants g^x+K- > 9 'f w+w- > 9 f K+K- can De estimated by using the information 
on corresponding decay rates B\ 
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9(f>K+K- 
47T 



1.66, 



9 f Q 1T+T\- 



0.105, 9 f° K+K ~ = 0.6 GeV 2 

47T 



(13) 



and I (a, b) is the function known in the analytic form PflJlTIl : 



I(a,b) 



2(a - b) (a - 6) 2 
The functions f(x) and ^(x) are given by: 



+ 



g{b- x )-g{a^) 



(14) 



/(*) 



arcsm 



2 1 



2V^ 



x > 



Z7T 



a/4x — 1 arcsin 
iyr^4^[ln2±-i7r] 



a; > 
x < 



i 

i' 
i 

I' 



(15) 



with r)± = (1 ± VI -4x)/(2x) 



5. C7wraZ [/mtan/ Approach (U\PT) @JTT[. 

In this case the decay — > 7r + 7r~7 occurs through a loop of charged kaons that subse- 
quently annihilate into 7r + 7r~7. The fo resonance is generated dynamically by unitarizing 
the one-loop amplitude. Using notations of Ref. the form factor in Eq.(|I]) reads: 



U X PT 



^ch 



G V M^ 



Qi 



,2 ' 



+ 



K 



m K m 



K , 



x/2 



Fv 
2 



Gy G K + K - 



(16) 



where the coupling Gy and F v are related to the decays — > K + K~ and — > e + e~, 
respectively, f n is the pion decay constant and I(a,b) is the function given in Eq. (|i~4l) . 
Gr+k~ is defined by the integral: 



( ]max 



2 * 2 J ° Jq 2 + m\ (Ql - 4(g 2 + m 2 K ) + ie) 



Gr+k- 



In Eq. ([T6|) t c h is the strong scattering amplitude 

1 



q 2 dq 



(17) 



ch 



1=0 

KK, 



(18) 



The scattering amplitude t^g- is determined by using chiral perturbation theory (see 
Ref. @). We have used the following values for the above constants: Gy = 0.055 GeV, 
F v = 0.165 GeV, f w = 0-093 GeV, g max = 0.9 GeV. 

In Fig.|l] we present a comparison of the photon spectrum obtained using the event 
generator and retaining only the term |.M^| 2 in the cross section (cf. Eq.(|9])), with the 
analytic expressions from Refs. |3],|5]||. One sees a good agreement between the Monte Carlo 
simulation and the analytic results. Note also, that different models predict different shapes 
of the photon spectrum. 
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III. STUDYING THE DIRECT DECAY 4> -> vr+vr 7 AT DA$NE 

We now address the question of whether precision studies of the direct decay — > 7r + 7r~7 
are possible at DA$NE. While writing the general formula for the process e + e~ — » 7r + 7r~7, 
we have pointed out that the observable signal of — > 7r + 7r~7 might strongly depend on the 
interference with the FSR. The fo signal may be enhanced if the sign between l-M^I 2 and 
2Re {.M^JW^} is the same (constructive interference) or may be reduced in the opposite 
case (destructive interference). 

In Fig.^| we present the spectrum of photons in the reaction e + e" — > 7r + 7r~7 in the 
situation when the invariant mass of two pions is close to the mass of the fo meson. We 
consider both constructive and destructive interference and generate events with and with- 
out collinear radiation, but initial and final state radiation is always kept. One sees from 
FigfJ that the collinear radiation results in the reduction of the signal. However, if the 
tagged photon is emitted at a relatively large angle, the effect of collinear radiation can be 
partially removed by combining the information on the position of the neutral cluster in 
the calorimeter with the directions of the charged pions determined with the drift chamber. 
In this case the kinematics of the reaction becomes over- const rained and it is possible to 
restore the "actual" center of mass energy for any given event. This will require a dedicated 
analysis, however. 

Fig.| shows the signal-to-background ratio 

S/B = (19) 

d(TQED 

for different models. As expected, the sign of the interference affects not only the magnitude 
of the decay 4> — > 7r + 7r~7, but also the shape of the distribution. The models where the 
structure of the fo meson is assumed show a broader signal for the constructive interference^] 
than in the opposite case. In addition, the "no structure" model does not show a clear 
peak in the case of destructive interference and the U\PT model in the case of constructive 
interference. 

The number of events required to separate the signal from the background can be ob- 
tained by estimating the necessary number of events in the energy region around the fo 
peak. We require the statistical error to be smaller than 10% of the signal itself. Hence, 

AN 

5N < — , (20) 
~ 10 ' V ' 

where AN is the number of events due to direct decay of the <fi meson: 

AN = Atotal — A^QED • (21) 

If we introduce a parameter £ such that 

iV t otal = (l + O^QED, (22) 



2 The excess of events is significant in the region of photon energies 20 MeV < E~ < 100 MeV. 
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TABLE I. The number of events and the integrated luminosity required to observe the direct 
decay (ft — ► 7r + 7r~7. 



Eq.(p0|) takes the form: 

5N < (23) 
~ 10 K ' 

The value of £ can be estimated from the S/B ratio shown in Fig.|3]. Using the standard 
formula for statistical fluctuation, 5N = \/N tata x, we estimate the number of events required 
to separate the contribution of (ft — ► 7r + 7r~7 from the QED background: 

AW>(1 + 0(T) 2 , ^>^- (24) 

Here e is the overall detector efficiency for 7r + 7r~7 events and C is the required integrated 
luminosity. 

The results are summarized in Table |. We use e = 0.5 and ctqed = 2.1 nb with the 
cuts I cos# 7 | < 0.9, | cos^±| < 0.9, | cos^l < 0.9 ( without collinear radiation). When the 
collinear radiation is included, ctqed is reduced to approximately 2 nb. 

In a similar way we make a rough estimate of the number of events and the luminosity 
required to discriminate between different models for the direct decay. We obtain: 
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TABLE II. The number of events and the integrated luminosity required to distinguish between 
two different models for the rare decay <j> — > ir + ir~'y. 



n qed > £ 



(25) 



with 



62 = 6-6, (26) 

and 6,2 are the ^-parameters for two models under consideration. In Table [TT] we summarize 
the results. One can see that the required number of events can be accumulated at DA$NE 
in less than one year assuming the luminosity £ = 10 31 cm -2 sec -1 . The required luminosity 
is, however, only indicative. 



IV. DIRECT DECAY OF THE cf> MESON AND THE MEASUREMENT OF THE 
ELECTRON POSITRON ANNIHILATION CROSS SECTION AT DA$NE 

It was suggested in Ref. H that the measurement of <r(e + e - — > hadrons) at DA$NE 
for different values of the center of mass energy can be performed by analyzing events with 
additional hard photon emitted at a relatively large angle (9 1 > 7°). The difficulties of this 
approach are related to the obvious fact that the hard photon can be emitted from both 
initial and final state of the process. If the ISR takes place, the total energy of the collision 
is reduced and such events can be used to measure a(e + e~ — > hadrons) at different energies. 
In contrast to that, the photons caused by the FSR represent a background that must be 
suppressed by applying suitable cuts. Since, to be competitive |],|2[], the measurement of 
a(e + e~ — > hadrons) for y/s < 1 GeV has to be performed at the one percent level, the 
practical realization of this idea is a non-trivial experimental task. In Ref. H only the QED 
process was studied. Here we would like to add the direct decay — > 7r + 7r~7 which also 
contributes to the FSR and therefore increases the background. In Fig.^ we show the values 
of (dcr to tai/d-E 7 )/(d(TisR/d-E , 7 ) with and without the contribution of the direct decay. The 
"pure QED" case was studied in Ref. @. The photon energies 20 MeV < E 1 < 100 MeV 
are considered, which corresponds to 0.836 GeV 2 < Q^+ w - < 0.996 GeV 2 . These invariant 
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masses of two pions include the contribution of the fo resonance and for this reason the 
largest contribution of the direct decay is expected in this region. The cuts reduce the FSR 
considerably; nevertheless, its contribution close to fo peak is significant. 

As discussed in 0, even "pure QED" theoretical predictions for the FSR are, strictly 
speaking, model dependent. It is therefore important to get a handle on it experimentally. 
In Ref. 0] it was suggested to use the forward-backward asymmetry of the produced pions 
to control the FSR. The direct decay <ft — ► 7r + 7r"7 changes the forward-backward asymmetry 
in the expected manner. Since the contribution of the direct decay is significant only if the 
invariant mass of the two pions is close to the mass of the fo meson, the forward-backward 
asymmetry integrated over large range of Q^+ n - is not affected by the direct decay. Hence 
it can be used to control the models for QED-like final state radiation. On the other 
hand, by applying the cut 0.836 GeV 2 < Ql +7T - < 0.996 GeV 2 , we significantly enhance 
the contribution of the direct decay to forward-backward asymmetry. This is shown in 
Fig.|5] where predictions of K + K~ model are displayed for both constructive and destructive 
interference. 

V. CONCLUSIONS 

We have discussed the contribution of the direct decay — > 7r + 7r"7 to the process 
e + e" — > 7r + 7r~7 at DA$NE energies. To facilitate this study, three different models^ for 
the direct decay <fi — > 7r + 7r"7 have been implemented into the Monte Carlo event generator 
described in Ref. @. 

The importance of this decay is twofold. First, it gives the information about the nature 
of the / (980) meson. Second, it provides an additional background to the measurement of 
a(e + e~ — > 7r + 7r") at different values of the center of mass energy by tagging the hard photon 
in the reaction e + e" — > 7r + 7r"7. 

We have shown that DA<3>NE has a very good potential to study the nature of fo reso- 
nance. Even with moderate luminosity £ = 10 31 cm" 2 sec" 1 , it is possible to discriminate 
between different models for the decay <p 7r + 7r"7 in a relatively short time. 

As for the measurement of the hadronic cross section a(e + e~ — > hadrons) at y/s < 
1 GeV using the process e + e~ — > 7r + 7r~7, we have found that the direct decay 4> — > 7r + 7r"7 
increases the final state radiation by several percent in the region of pion invariant masses 
0.836 GeV 2 < Q 2 +^_ < 0.996 GeV 2 , but quickly dies out beyond this region. 

Finally, we note that it will also be possible to perform a detailed study of the decay 
(p — ► 7r°7r°7 at DA$NE . We believe that it will be quite useful to combine these independent 
measurements in order to check the theoretical understanding of the fo meson. 



3 It is relatively straightforward to include other models for the direct decay, for example the four 
quark model of Ref. ||, to the event generator. We plan to do that in the nearest future. 
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FIG. 1. The energy spectrum of photons produced in the direct decay eft — > ir + ir~^f without 
initial and final state radiation. The results of the Monte Carlo simulation are compared with the 
results of analytical calculations (solid curves) for different models |]|^||. We use rrif = 0.976 GeV, 
Tf =34 MeV. For the "no structure" model the branching ratio Br(^> — > /07) = 10 -4 has been 
used as an input (this branching ratio is of the same order of magnitude as is predicted by the 
other two models). For the purpose of comparison, no exponential damping has been applied to 
the "no structure" model and the events were generated without collinear radiation. We have 
applied the following cuts on the polar angle and the energy of the photons: 10° < # 7 < 170°, 
20 MeV < £L < 100 MeV. 
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FIG. 2. The spectrum of photons in the reaction e + e — > 7r + 7r 7. The following cuts were 
applied: |cos# 7 | < 0.9, |cos0„.±| < 0.9, | cos^^l <C 0.9, where 6^ and 9 n are the polar angles of 
photon and pions, respectively, and 6^ is the angle between the photon and the 7r + momenta 
in the center of mass frame of the two pions. The exponential damping was applied to the 'no 
structure'" model and the branching ratio Br(0 — > /07) = 2.5 x 10 -4 has been chosen. Different 
models are distinguished by different hatching. "Pure QED" means that only the contribution due 
to dciQED is considered. 
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FIG. 3. The signal-to-background ratio (dctotai/d-E 7 )/(d<TQED/d.E 7 ) as a function of photon 
energy. The cuts are | cos6> 7 | < 0.9, | cos^ij < 0.9, | cos 6^1 < 0.9. See text for more details. 
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FIG. 4. Ratio (d<7 to tai/di^ 7 )/(d<7isR/d.E 7 ) as a function of the energy of the photon. The 
cuts are 7° < # 7 < 20°, 30° < 6 n < 150°, and the invariant mass of detected particles in the 
final state Q 2 +7r - 7 > 0.9 GeV 2 . "Pure QED" ratio is defined as (daisR+FSR/di? 7 )/(d<7isR/d£' 7 .) 
Different pictures correspond to different models for the fo resonance. The cases of constructive 
and destructive interference are considered. See text for more details. 
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FIG. 5. 7r + angular distribution. The photon angle is 60° < # 7 < 120°. The invariant mass of 
the two pions is 0.836 GeV 2 < Q 2 +7r _ < .996 GeV 2 . 
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